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CALCULATION OF LINEARIZED SUPERSONIC FLOW OVER 
SLENDER CONES OF ARBITRARY CROSS SECTION 

By Vincent R. Mascitti 
Langley Research Center 

SUMMARY 

Supersonic linearized conical -flow theory is used to determine the flow over slender 
pointed cones having horizontal and vertical planes of symmetry. The geometry of the 
cone cross sections and surface velocities are expanded in Fourier series. The sym- 
metry condition permits the uncoupling of lifting and nonlifting solutions. The present 
method reduces to Ward’s theory for flow over a cone of elliptic cross section. Results 
are also presented for other shapes. Results by this method diverge for cross-sectional 
shapes where the maximum thickness is large compared with the minimum thickness. 
However, even for these slender-body shapes, lower order solutions are good approxi- 
mations to the complete solution. 


INTRODUCTION 

The solution to supersonic flow over a cone of arbitrary cross section has been 
treated extensively in the literature. References 1,2, and 3, for example, treat this 
problem with nearly exact inviscid formulations, which admit the existence of shocks and 
vortical singularities induced by large crossflows. The present solution to linearized 
supersonic flow over a cone of arbitrary cross section can be superimposed to obtain the 
flow over a body which changes shape longitudinally. An example of this approach is the 
classical solution for the body of revolution obtained by the superposition of circular - 
cone solutions. (See ref. 4.) The superposition approach could possibly be applied to 
the solution of the flow over wings, fuselages, and wing -body combinations. 

A general theory for the solution of slender bodies as well as cones was given by 
Ward in reference 5. Ward’s theory indicates that for slender bodies, the velocity poten- 
tial satisfies Laplace's equation in the cross -plane coordinates; and thus, the methods of 
classical hydrodynamics can be used to obtain solutions. Ward's theory has been used in 
detail only to obtain the flow over an elliptic body (refs. 6 and 7), since the transformation 
from ellipse to circle by the Joukowski transformation is well known. 

The purpose of this paper is to present a method for determining the flow over a 
cone of arbitrary cross section based on the assumptions of supersonic linearized theory. 


The method of solution will not depend upon incompressible cross flow (slender -body 
theory); therefore, results by complete linearized theory can be obtained. The method 
uses the conical -flow solutions to the wave equation first classified in reference 8. 

These solutions individually have the property of vanishing at the Mach cone. The cone 
geometry and velocities on the surface are expanded in a Fourier series to satisfy the 
surface boundary condition. 

The validity of this method is demonstrated by comparing the present results with 
those from Ward's theory for the elliptic cone. Results from the present method are 
shown to diverge when the maximum thickness is large compared with the minimum thick- 
ness. Results are also presented for two cones of arbitrary cross section with slender- 
body approximations. Finally, the present method is applied to the elliptic cone without 
slender -body approximations, and results are compared with an extension of Ward's 
theory presented in reference 9. 

The method of this paper is presently restricted to cones with horizontal as well 
as vertical planes of symmetry. With this restriction, the lifting and nonlifting solutions 
can be treated separately and superimposed to obtain a complete solution. 

SYMBOLS 

Ajj kth coefficient of body-geometry expansion 

a,b cone cross-section parameters (see fig. 1) 

B(x,r,0) function describing body surface 

Cn nth coefficient in a series of superimposed potential solutions 


G 


M 

m number of terms required to approximate body geometry 

N total number of superimposed solutions 

n individual solution to be superimposed 


pressure coefficient 


k=°° 


body -geometry function, G = cos (k0) 

k=0,l,2 

free -stream Mach number 


2 



kth coefficient of boundary equation 


velocity vector in flow field 


kth coefficient of expansion, cos nd G n 
kth coefficient of expansion, sin nd G n 


free -stream velocity 


body -axis system 


x,r,0 


cylindrical coordinate system (see fig. 1) 


angle of attack, deg 


|8 = V m<S ‘ 1 


perturbation velocity potential 


Subscripts: 


k,m,n integer indices 


METHOD OF SOLUTION 


Governing Equations 

A uniform supersonic stream of Mach number M flowing over a cone at an angle 
of attack a with respect to the X-axis (fig. 1) is considered. Under the restrictions of 
inviscid flow and small velocity perturbations, the flow field can be described in terms of 
the perturbation velocity potential <p which obeys the equation 

8r 2 r 8r r 2 802 P ax 2 

The general solution for the flow over a pointed body is obtained from reference 10, 
equation (18), as 


0(x,r,0) = - ^ 


cos nd 


n=0,l,2 


'P r F n (£)cosh(n cosh -1 ^4)d£ 


l(x - I) 2 - /3 2 r 2 


3 


where F n (£) are functions to be determined by the surface boundary condition. If 
4 = x - j3r cosh z, then d£ = 1 - /3r sinh z dz and since properties along a ray are con- 
stant for conical flow, 

F n (£) = F n £ = F n (x - /3r cosh z) 

Then the general solution for a cone becomes 

n=°° 

* (x ' r - e) - h 2Z. 

n=0,l,2 


r° 

F n cos n0 \ . ( x - /3r cosh z)cosh nz dz 

J , -1 x 


'cosh- 1 ^ 


(1) 


The integrated forms of equation (1) were first given individually in reference 8 as 
solutions to conical flow having the property that the perturbation potential vanishes at 
the Mach cone. These perturbation velocity components are 


n=°° 

&..1- > F n 


f 


0 


cos n0 \ „ cosh nz dz 

* 2 %%72 “ 


(2a) 


8r 


11=00 0 

-i_ > F n cos n0 \ cosh z cosh nz dz 

2tt 4— _ x 


(2b) 


n=0,l,2 


and 


n=°° 


1 96 _ 
r 80 


cosh_1 ^ 


-1— y F n n sin n0 \ . (x - 0r cosh z)cosh nz dz 

2nr 4— „ n 


n=0,l,2 


'cosh" x ^ 


(2c) 


Integrating equations (2) yields 

n=°o 

^ = -CqZ - cos 0 sinh z - C n cos n0 

n=2,3,4 


(3a) 


-2^ = C o 0 sinh z + Cj cos 0 [jcosh z sinh z) + zj + ^ 


n=2,3,4 


rn- n( 4~ sinh ( n + ^ Z i sinh < n ’ ^ 

COS ne L _ 2(nTI ) + 2(n - 1) 


(3b) 


and 


r = 2 C 1 Sin *G C0Sh z sinl1 z ) - Z J - 0 


n=°o 


nC n sin n0 


n=2,3,4 


sinh(n + l)z sinh(n - l)z _ cosh z sinh nz ] 
2(n + 1) 2(n - 1) n J 


(3c) 
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Ill 


where 


and 



z = cosh 


-1 x_ 

0r 


With the appropriate boundary condition, the term n = 0 is the required linearized 
theory solution for the nonlifting circular cone and the term n = 1 is the solution for the 
lifting circular cone (ref. 4, pp. 215 and 223). The infinite series represents a general 
solution to a cone with an arbitrary cross section. It is expected that as the cross- 
sectional shape approaches a circle, a lesser number of solutions will be required. 


If the cone is slender, that is, 

x 
fir 


£-» 1 


then 


z = cosh’ 1 * In £ 


P r 


pr 


and 


sinhz=l-fl--l*f- 
V/3 2 r2 Pr 


cosh z = — 

pr 


With the aid of hypergeometric identities, equations (3) become 

n=°o 


|£ = _c 0 Inf* - 
9x u /3r 


n + 1 


n=l,2,3 


n C n cos n^- 


n 


(4a) 


n=°° 


1 * = C o(f) + C 1 cos 6 (f) + ^ c n cos n ®(|) 

n=2,3,4 


n+1 


(4b) 


and 


n=°° 


lie = C 1 sin e (?) + S Cn sin ne (!) 

n=2,3,4 


n +1 


(4c) 
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In agreement with the slender -body theory of Ward (ref. 5), the velocity components 
and i satisfy Laplace's equation: 


3 r 2 r 96 r 2 m 2 


Surface Boundary Condition 

Since any body shape can be described as B(x,r,0) = 0, the surface boundary con- 
dition for steady flow is given by q • VB = 0 where q is the velocity vector. Com- 
bining these expressions gives 

(u cos a + + fu sin a cos 0 + 

\ ax /ax V ar / ar 

+ [-U sin a sin 0 + i = 0 

\ r 80/ r 80 

where all quantities are evaluated at the body surface. For flow over a cone 

B = r - — = 0 
G 

and the boundary condition becomes with rearranging 
U cos a + = (u sin a cos 0 + &J G 


With the assumptions that a is small and U cos a » and with the perturbation 
velocities normalized by U, 


l-«fcos0G-sin0§)=G^ + |^I|^ 

\ 80 / 3r 30 r 80 


If G is expressed by 

k=°° 
G 


= * = A k cos k0 

k=0,l,2 


(5) 


the velocity components (eqs. (4)) can be expanded in Fourier series and the boundary 
condition can be put in the form 
k=°° 

Pjj cos k0 = 0 

k=0,l,2 
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Since cos 06, cos 9, cos 29, . . cos k0 are linearly independent, setting Pq = 0, 
Pj = 0, . . ., P k = 0 will give k linear equations for the constants Cq, Cj_, 

• • •> Cfc. 


If 


G = 


k=0 

y~ A k 


k=0,2,4 


cos k0 


( 6 ) 


then the cone cross section is horizontally as well as vertically symmetric. The n-even 
velocity contributions contain only even cosine terms, and the n-odd velocity contribu- 
tions contain only odd cosine terms. Rewriting the boundary condition (eq. (5)) gives 


- «(cos 9 G - sin D §) = (o * (§ ) If) ♦ (° !f) oa(i * (f ? 

V . — .✓ V — — . ^ — — . — — — 


Constant 


Odd cosine series 


Even cosine series 


Odd cosine series 


Therefore, with the assumption of equation (6), the n-even solutions are nonlifting solu- 
tions and the n-odd solutions are lifting solutions, and the two problems can be solved 
separately and independently. The nonlifting solutions obtained from equations (4) are 

n=°o 


X ^CnCOsneG" 


n=2,4,6 


n=oo 


= / C n cos n0 G n+1 

9r n 


(7a) 


(7b) 


and 


n=0,2,4 


n=oo 


? C n sin nd G n+1 
n=0,2,4 

where the boundary condition from equation (5) is 
dr 90 r 90 

The lifting solutions obtained from equations (4) are 
n=°° 


(7c) 


(7d) 


90 _ n + 

ax n 


- cos n6 G n 


(8a) 


n=l,3,5 
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(8b) 


and 


n=°° 


J£= C n cos n0 G n+1 


n — 1 j3 j5 


n=°° 



n=l,3,5 


C n sin n0 G n+1 


where the boundary condition from equation (5) is 


(8c) 


-a(^ 


cos 0 G - sin 6 


88 j 9r 90 r 90 


The ability to split the boundary condition into two parts also holds for equations (3) 
and has previously been demonstrated only for the circular cone in reference 11, page 241. 
With slender -body assumptions, the nonlifting solution contains a Mach number variation 
only in the first term of the 90/ 9x expression, which does not enter into the solution of 
the boundary condition. The lifting solution is entirely independent of Mach number. 


Expansion Procedure 


There remains the necessity of stating the procedure by which the infinite series 
are truncated and a finite solution obtained. Results of the expansion must approach 
those obtained by Ward’s theory as n approaches infinity. Since the expansion proce- 
dures for the lifting and nonlifting solutions are similar, only the procedure for the non- 
lifting solution will be given. 


k=m 


(1) Expand G in the form 


Ak cos k0 


k=0,2,4 

(2) Choose n = 0, n = 2, n = 4, ..., n = N number of terms and expand the 
velocity equations (eqs. (7a) to (7c)) in the form 


n=«> 


- 

dx 


r in 

-Co 1“ — 


n=N 


k=N 


n + 1 
n 


-n 


n=2,4,6 


cosnOG"--]^ ^C n > 


n=0,2,4 


k=0,2,4 


R k , n cos k0 


n=°° 


n=N 


k=N 


2Z CnCOsneG' 1 ^ ^ C n R M 


cos k0 


n=0,2,4 


n=0,2,4 k=0,2,4 
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and 


n=°° ji=N n=N 

g?|£ = C n sin n0 G n = C n 2Z S k,n sin k0 

n=0,2,4 n=0,2,4 k=0,2,4 

Although the function cos n0 G n is completely defined by a series to n(m + 1), it is 
truncated at N. 

(3) Solve for the boundary condition in the form 

1 = G 2 (h^-] + i — (g 2 ) ( 

1 ^ \U dr J + 7 ae K ' \Gr aej 

(4) Solve for pressure coefficients by 



This expansion procedure is best illustrated by presenting, in detail, the solutions 
for N = 0 and N = 2. These solutions are presented in the section entitled "Nonlifting 
Solutions." 


RESULTS AND DISCUSSION 

The expansion procedures described previously must converge to well-known solu- 
tions based on the same assumptions. Therefore, a comparison of results from this 
theoretical method for the elliptic cone with those from well-known solutions is an impor- 
tant test. 


Nonlifting Solutions 

Elliptic cone .- As an illustration of the expansion procedure, the lowest order solu- 
tion (N = 0) and the next to lowest order solution (N = 2) will be presented in detail for the 
elliptic cope. The geometry function G for the elliptic cone is 


G = 7 = V A 0 + A 2 cos 26 

where 


a 2 + b 
2a 2 b 2 


A 2 = 


b 2 - a 2 
2a 2 b 2 


and a and b are the semimajor and semiminor axis, respectively. For N = 0 the 
velocity components obtained from equations (7a) to (7c) are 


= -C 0 In = -C 0 In 
9x 0 j3 0 


a + b 
j3ab 
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10 


( 9 ) 


this equation becomes 


( 10 ) 



or 


4aV 
P ' a 2 + b 2 


Cn - 


In 


a + b 
/3ab 


-1 + 


2a 2 b 2 


1 + 


a 2 - b 2 


a 2 + b 2 V a 2 + b 2 


cos 26) 


For a = b (slender circular cone), 


C n = 2a 2 1 In 


/3a 2 


(ID 


( 12 ) 


Equation (12) agrees with that presented in reference 11, page 234. Ward’s slender-body 
theory as given by Van Dyke in reference 9 is 


Cp = ab 


2 In 


/3(a + b) 


2 + 


ab 


a 2 sin 2 ?7 + b 2 cos 2 ?7 


(13) 


where 


tan rj - — tan 0 
b 


This expression can be shown to have the following Fourier series expansion in 0 (see 
ref. 9): 

k=°° 


Cp = 2ab 


In 


4 


- 1 


0(a + b) 

For k = 2, the expression becomes 

4 


2a 2 b 2 'X /a 2 - 

a 2 + b 2 ,/n o"Aa 2 + 


k=s0,2,4' 


-bf 

b 2 > 


Icos k0 


C p = 2ab 


In 


/3(a + b) 


- 1 


2a 2 b 2 A a 2 - b 2 
a 2 + b 2 \ a 2 + b 2 


cos 20 


(14) 


The expression for the present method (eq. (11)), where N = 0, is repeated here for 
convenience: 


Cp 


4a 2 b 2 
a 2 + b 2 


In 


a + b 
/ 3ab 


-1 + 


2a 2 b 2 
a 2 + b 2 


1 + 


a 2 -b 2 
a 2 + b 2 


cos 20 


Comparison of the present method (eq. (11)) with Van Dyke's solution (eqs. (13) and (14)) 
is shown in figure 2 for M = \/2, x = 1, a = tan 30° = 0.57735, and b = 0.57735, 0.5, 0.3, 
and 0.1. Although the geometry shown in this figure is far beyond the range of linear 
theory, and certainly slender-body theory, the purpose of this figure is simply to compare 
theoretical results. 

In the present method for N = 2, the velocity components (given previously as 
eqs. (7)) are 


11 


80 

8x 


-Co In - | C 2 COS 20 G 2 

- c o( ln isr - rrl cos 2e ) - 1 c 2 (t + A » cos 29 ) 


(15a) 


and 


!!£= C n + Co cos 20 G 2 
G 8r u z 

= C 0 + C 2 ^ + A 0 cos 20 


-L = o + Co sin 20 G 2 
Gr 80 2 


= 0 + C 2 Aq sin 20 

Substituting equations (15) into the boundary condition (eq. (7d)) yields 

1 = (A 0 + A 2 cos 20) (c 0 + + C 2 A 0 cos 20^ 

- A 2 sin 20(C 2 Aq sin 20) 

By expanding equation (16) 

^2 a 2 


1 " AqiCq + 


2 . * . C 2 A 2 


~ \ c 2 A 0 + a 2 c 0 + 


cos 26 


- I C 2 + C 2 A 0 A 2 jcos 40 = 0 

is obtained, and the boundary condition is now in the form: 
k=N 

^ Pjj cos k0 = 0 


k=0,l,2 


Therefore, 


Co An 

Pq = ° = 1 - a oI c o + ~t 


(15b) 


(15c) 


(16) 


p 2 = o = c 2 a 0 2 + a 2 c 0 + c 2 ~ y ~ 
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or 

Ao 

and C 9 = - — =- 

2 An 3 
A 0 

A comparison of the present method for N = 2 with Van Dyke's solution is presented in 
figure 3. 

The expansion procedure just illustrated has been generalized and programed so 
that higher order solutions can be obtained. Figures 4 and 5 show results for N = 8 
and N = 32, respectively. Although the calculation procedure does not converge for small 
eccentricity, lower order solutions are good approximations to the complete solution. 

Arbitrary cone .- The present method has been used to calculate the pressure dis- 
tribution over cones with the geometry shown in figure 6. These results are shown in 
figures 7 and 8 for N = 8 and N = 32, respectively. The results for the most winglike 
cone are diverging. 

Lifting Solutions 

The lifting solution has been shown to be independent and separable from the non- 
lifting solution. Solutions are obtained by superimposing n-odd solutions in the identical 
manner as was done with the nonlifting solution. Again the agreement between the results 
from the analytic solution to the lifting elliptic cone and those from the present method 
will be the important test of the expansion procedure. It should be noted that the lifting 
solutions are independent of Mach number. 

Results for the elliptic cones previously examined are presented in figures 9 and 10 
for N = 7 and N = 31 at a = 10°, respectively. These results also show signs of 
divergence in the thin-wing limit. 

Figure 11 shows results for the cone of arbitrary cross section at a = 10° for 
N = 32. Large changes in pressure are shown in the region where the straight-line geome- 
try changes to circular geometry. This trend is also indicated by the experimental data 
of reference 12. 

A listing of the computer program to calculate the pressure distribution around lift- 
ing cones of arbitrary cross section with slender -body theory is presented and discussed 
in the appendix. Computational time for obtaining pressures with this program has been 
estimated at 1 minute per case for N = 32 on the Control Data series 6600 computer 
system at the Langley Research Center. 


c o = fl 

A 0 
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Results Without Slender -Body Assumptions 

The present method of Fourier series expansion can be applied to the solution of the 
cone with arbitrary cross section without slender-body assumptions. The velocity com- 
ponents, in this case, satisfy equations (3). Results from this solution are presented in 
figure 12 for several elliptic cones. The solid line represents results obtained using the 
"not-so-slender" solution of reference 9. The solution of reference 9 does not reduce to 
the well-known circular -cone solution of reference 4, page 214. However, the present 
method does reduce to identically the correct circular -cone solution. 

Theoretical Limitation 

For all the results shown, the present method diverges in the thin -wing limit; how- 
ever, the lower order solutions gave good approximations to the actual solution. 

Ward, in reference 5, indicates that slender-body theory should not be applied to 
cross-sectional shapes where the local radius of curvature is small compared with the 
maximum thickness. This restriction could explain the divergence of the present method 
as the geometry approaches a wing. It is important to recognize that the analytic solution 
to the elliptic cone given in reference 6 was achieved only after the transformation from 
ellipse to circle by the Joukowski transformation. The results of the present paper would 
indicate that in order to achieve converged results for arbitrary winglike cross sections, 
an initial transformation would be necessary. 

CONCLUDING REMARKS 

A Fourier series expansion procedure has been developed to solve for the flow 
around slender cones at supersonic speeds. The results for an elliptic cone reduce to 
those obtained by Ward's theory. Both lifting and nonlifting solutions to the arbitrarily 
shaped cone can be obtained, except in the thin-wing limit where the method diverges. 

The present method has been programed and applied to cone solutions without slender- 
body assumptions with good results. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 12, 1972. 
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APPENDIX 


COMPUTER PROGRAM FOR CALCULATING THE PRESSURE DISTRIBUTION 

AROUND SLENDER CONES 


The calculation procedure described in the main body of the paper for obtaining 
pressure distributions around slender cones has been programed for high-speed digital 
computation. The purpose of this appendix is to provide a description of the necessary 
input and available output as well as a FORTRAN IV (ref. 13) listing on the source pro- 
gram. An example input case and the resulting output listing are included. 


Description of Program 

The program reads in the necessary number of lifting and nonlifting solutions to be 
superimposed. The coefficients of the required geometry function are computed by 
standard Fourier expansion series techniques (numerical integration). The matrices 
given by the boundary condition are formed and inverted using Gaussian elimination. The 
values for Cn that result are used to compute the velocity components and pressure 
distribution. The program listing that follows has the geometry of an elliptic cone built 
in. However, this geometry can readily be replaced by any arbitrary description of 
geometry. 


Program Listing 

The FORTRAN IV listing of the source program used on the Control Data series 6600 
computer system at the Langley Research Center is as follows: 

COMMON THETA (l Pi > *0(181) * 1 1 * P I , Q0 . QQ ( 840 ) *A(17«18)*R(17)*K 
DIMENSION MlPl ) *RR0(32> * RO ( 32 ) * R ( 32 0 * 32 ) * 

1P(32), RR(3?*3i)* TT (32*32) iF (181) « H ( 18 1 ) * 

2FF ( 181 > *Hh < 181 ) *C (32) *CX (32) *0° (32) *CT (32) *CP < 1 Hi ) 

3* AR (32*32) *AS(32* 32) *AT (32*32) 
real m*lam 

NAMELTST/NUM1 /mm«Im»NW 
NAMELIST/NUM2/M, ALPHA 
READ (5*NUM1 > 

WRITE (6.NUM1 ) 

DO 10 I = 1 « N 
RO ( I ) =0 . 

LL=MM»N 
DO 10 J=1,LL 
10 R ( J* I ) =0 • 

P 1=3* 1 41 592.653589793 
AA=. 57735 
66= . 3 

theta <n=n. 

DO 1 1=2*181 
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APPENDIX - Continued 


1 THFTA ( I ) =THET A { I — 1 > ♦ t . 

DO 2 1=1,181 

thetr=theta <n i/iro. 

2 G ( I ) =SORT ( ( AA**2 + BB<H>2* <BH**2-a A»* a) *C0S <2.*THFrR> ) / <2.*AA**2*8B** 
12) ) 

100 FORMAT (2X1F16. 8) 

DO 3 K=l ,? 

DO A- I = 1 * T 8 1 

4 Q ( D =o ( i ) 

1I=MM*K 
CALL INTER 
HO (K) =00 

DO 5 J=2* 1 1 *2 

5 R ( J * K ) =00 ( J) 

3 CONTINUE 

DO 7 T=l,|8l 
7 0<I)=AL0G<2.*Gm ) 

1 1 sN 

CALL INTER 
R0=00 

DO 8 J=2»TI*2 


a i- ( J) -TO ( J ) 

',0 ?• L = I • N 
RW> (i 1 =0 • 

DO 2 * K = 1 • n 
?4 RR <K.( > = v . 

DO 2 I I = 1 4 1 a 1 
Tt-tETR=ThL T A < I >*RT/1*0. 

21 n ( I > =r05j <L*THhTR> *0 ( I ) **|. 

T1=N 

CALL INTER 
wRfj (| > sO,- 
i'O 2 > J=? » i t . 2 
?■? RR ( J. L) =NU ( n 
20 < ONT i DiJF 
no 4 • L= I 4 N 

00 4 A K = 1 4 i'l 

4a TT(K.|.)=a. 

DO 41 1 = 1 4 1 a] 

ThFTR = THL T A f I > *P| / 18.-). 

41 n (I ) rSIN (L*ThETR> *(i ( H **l 

t i=n 

CALL ONTi-R 
.10 4? J = r 4 I T .2 

42 TT1J.I )=>)0(.l) 

4 fi OONTtniJF 

>0 5 a I = 1 • 1 K 1 
THt T r - = THF TA(I)»R 1/180. 

M ( I ) i . 

F(I)-R0«2) 

K K = 2 r- , .»'4 

no 5 - k=i,kf 

H ( I ) -M ( I ) -K*R <K . 2) *SIN (K*THFTR) 
SO F (1 )-F (1) ♦PIK»?) »COS(k*T'JFT.R) 

S4 CONTINUE 

no Si 1 = 1 4 HI 
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THETPsTHi- 7 fl < I 1 * u t/lHf). 

FF ( I ) = P0 ( 1 » 

no 51 K a t » N 

6) FF ( I > =FF ( ]) *P (K. 1 > *C05 (K»THFTR> 

no Bp I = i * 1 -p 1 

52 <)(T) rf ( I > 

1I=N 

CALL INTER 
A ( 1 . 1 t =<K 

no 5n K=2 .it.2 
L. =K/p* 1 

53 A (L » T > =UU ( X ) 

00 61 J=2*N.? 

..IJ=J/p + 1 

no 6i i=i *i5i 

THET D = THF T A (I >*PI/1H0. 

F F ( T > = Rku (.J) 

HH(I I =0. 
no 61 K* 1 * N 

FF(Ii=FF(I>+RR(K«J)*CnS(K*THFTR) 
61 HH ( I \ rHH t I ) *TT (K * J) *SlN (K*THFTR> 

no 6 p i=l .] pi 

6 ? 0 ( 1 ) :F F ( I) *F ( I ) +. 6 *H (T )<MH( I ) 

1 I =N 

TALL 1 NT! R 
All, ).l) =Oo 
DO 6i « = ,-■, IT, 2 
L=K/p* 1 

63 A (L . , II) = uO ( « ) 

I = N/p'*2 

A ( 1 .1 ) =1 . 

no 6s J=2.N.? 

A=J/P+1 

65 A (K,| )=U. 

I=?«'*f.' + ? 

TF<I.GT.N>00 TO 67 
no 6* J= 1 » N • 2 

K = J/P ♦ 1 

66 A (K , 1 \ =0. 

no 6^ l=c*n .2 
I.L = L/p+1 
I =?*mm+MM»l+L+2 
TFd.GT.MRO TO 67 
no 6« J=I *N.2 

K-rrJ/p* 1 
6P A (K,|.( ) sj, 

67 CONTINUE 
k-=N/p + 1 
k = N/i* 1 

tall '^AT p 1 X 
CO = R M ) 
nO 7 ‘ L=?»N.2 
J=L/p+1 
70 0 ( L ) ■= R ( J ) 
no I?'. L= 1 «mW 
no 1 ?a K= 1 »MN 
12A ar(k.I )=j. 

no 1 pi 1=1.131 

THETw = THt-.T A ( H *PT/1 «0. 

121 0(I)-=nOS(L*ThfcTR)*G(T>**l. 

I T =N M 

TALL A1MTEP 
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no 1 •>? . 1 = 1 . mn . 2 

1?? AR(J.L)=00(.|) 

1?0 rONTtwUF 

no n:, L=1 »mM 
nO 1 K = 1 , *>:('■ 

1.34 A S ( K . I 1 =i . 

nO lnl 1=1. 1*1 

THt'T r, = T Hr T A (I >*01/180. 

1 31 0(1) =003 ( L*ThETO) *G ( I » **l 

II =N, 

CALL A INTER 
r ’0 In? J=i,smI>i,R 
13 ? < J.l ) =‘JU ( i) 

130 CONTtimUE 

MO 1 A L=1.mn 
no l uu K = 1 , mm 
\uu AT(K,|. )=u. 

no 1 a i T = 1 • l « 1 

THH.' = ThhTA)T)*?I/1kii. 

]4l 0( I ) r-STN I L*THETR) *0 ( T ) **l 

T T =N’ 

TALL A ON TER 
00 1a/ J=1,mm«? 

14? AT<J.L>= ><J< I) 

140 OONTiM'F 

no 1 M j=i.mm*? 

JJ= ( i + 1 ) / ? 

rO 1 ^ 1 1=1.) Ml 

THE T TH 1 - TAf I>*PI/180. 

FT ( T )=o. 

4 H(Hr!|, 
fO l,.l K= 1 • MI-1 

FF( I ) =FF ( I ) +AS(K, J) *cns (K*THFTR) 

1 bl HH ( 1 1 =HH ( I 1 + AT (K. J) *STN (X*ThFTW> 

00 If? 1=1, )8l 

16? 0(1) -FF ( I ) *r ( I ) + . 5*H ( T 1 *HH ( I ) 

TT=Nm 

CALL AINTER 
no If? K=1 ,nn*? 

L = ( K + 1 1/0 
163 A (L, |J) =UQ (K) 

SOQ READ (F,NUM2) 

T F ( EOF , 5 1 
Q WRITr (6.NUM?) 

L AM = c;QRT (M**?-l . ) 

PP = P^, + AU'b ( 1 ./LAM) 

XC0Nf=-1 ,*co*pp 

YCONf=CO 

nil Ra J = ?,M.? 

xCONc=-l .*F|.OAT (J+l ) /FLOAT ( J)*RRO ( J)*C(J> + XCOnF 
80 YC0Mr=C(01*PR0(J)+YC0ME 

TOONrrO . 

no 8? L=?,N.? 

OX <L ) =-l . *On*H (L) 
r P ( L i = 0 . 
rT(L)=0. 
no R? K=?,N.? 

ox < L > =-l .*F| OAT (K+ 1 ) /FLOAT (K> *R p (L.K > *C («) *C* (L) 
OR < L > =C (6) »pr (L.K> +CR (L) 
ft? FT (Li=C( K )*TT ( L , K ) +CT (L> 

00 9" LLL = 1 .181 
THET^THtTA(LLLl* p T/l J 0. 
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i)=XC r '\ l IT 
V=YC'1 nE*0<LL L) 

V = 0 • 

no 91 L = 2 tNt? 

ti=lUcx (L> *COS <L*ThETp» 

\j=\J + ru (L) <>CnS (L*ThfTR> *0 (ILL) 

91 VarW + rT(L> *S T N ( L* ThF TR ) *0 ( LLL 1 
rP(Li L>=-2.*(U+(V**2+.j**9)/?.> 

90 rONTrMUE 

alPh°=ALPH A*PI / 1 HO • 
no l'l I = l*ml 

THETc=TH£TA ( I > *oi/le»0. 

101 n ( I > - r- ( I ) *Cns (THETP) - ,S»nTM(THFlV) *h ( I ) /o ( T ) 
r i=nm 

CALL AlNTER 
|. = (Nm+ 1 >/2+l 

no 1 t ? j= 1 .mm«? 

K= < J* 1 ) /2 

10? A (K ,| ) 5=-l ,*aLPHR*00 * J) 

K= ( Nm* 1 )/c 
CALL MATRIX 
oO l-xo L=I*mm*? 

.1= <L + 1 ) /d 
170 C(L>=m(J) 

nO 1 = p L = 1»N>N»2 
rx(L>=-?.*AJ(L.l)*CU) 

CP <L> =C ( 1 > *AS (L * 1 ) 

CT (Li =C( n *AT (L* 1 ) 
no 1°? K=3*mM«2 

CX(L>=-1.*(*+1>/«*AP<|_»K»*C(> < '>+CX(L> 

C« <L) =C <" > *AS (L**) +C9 (L> 

18? CT(U=C(K)*AT(L«K)+CT<L) 

'•19 IT" (6»20?> 

?0? FOHMaT ^X^HTHETA. 10/6HCP. A=0 .9X1 1 HCP« A-C.P . A = . : . 7 X 7HCPTOT A L / ) 
no 1 Q 0 LLL = 1 .181 
T HE To = T HE F A (LLL> *PI/1 «0 . 
aU=0. 
an/ = 0 . 

Awl=0 . 

no 1 qi L= 1 . mm *2 

Al)=Atl*CX (L) *C0S(L*THET9> 

av = A'/ + CR(L> *C0S(L*THET9) *R (LLL) 

191 A'V = Aw + CT(L)*SIN(L*THET9)*G(LLL) 

CPA=Al PH9**9-2.*AU- ( ALPH9*C0 S (ThETR) +A V) **5- ( ALPHW*STM ThFTP) 
1 — A W ) 0*2 
CPT=CP A + CP ( I LL) 

190 wPITr (6.201 ) THET A (LLL) .C°(LLL) .CPA.C-PT 
?01 FQRMaT (?X<tFif,i8' 

00 To 500 
?00 STOP 
FMO 
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SUBROUTINE INTER 

COMMON THET 0(181 I ♦0(I«l) . 1 1 * P I * QO . 00 ( 640 ) .Ml 7. H) <H I 17) .K 
OIMEMSION F ( l 81 ) 

areb=o. 

00 I 1=2*181 

1 AREB = AREH*-(0(I>+Q(I-lt )/?. 

00=APFB/180. 

no 2 J=2 » 1 1 * 2 

AREB= 0 « 

F(l)=Om 
00 3 1=2*181 

thetp=theta ar*Pi/i«o. 

F(I>-0(I)*C0S<J*TuFTWi 
3 ARER=flREB+(F (I) +F ( I~1 ) 1/8. 

2 00( J>=?.*ARFH/lflO. 

RFTURN 

END 


SUBROUTINE oNTER 

COMMor. T w £Ta. ( 1 8 It .0 < 1 >*1 > • I I ^ I "NO *'Ji> < 6-*0 » * a (i 7 . j 4) . h f m . f 
~l I ME"'S I ON FI]H|) 

.'■>0 2 1=2*11.? 

APER = I‘>. 

F(l)— 0(1) 

1)0 3 1=2*781 

THE T -> = THi: T A ( T ) *P 1/180, 

E (I)=0(I)*STN(J*THETR1 

3 ARERrrflRE H+ (F(I) + F (I — 1 > )/ J . 

? 00<J*=?.*ARFH/I«0. 
b F T U h 'i 

FimO 

SUBRollTIUF A T R I X 

roNMfvi THE T a { 1 8 l • • U • l a i ) . 1 1 • p l • u j « Uu’ ( E>*+ , J i *4 ( l 7* (nt *ni j / )»•*' 
l.=K*l 

o.O 4 T T = 1 *K 
-•0 1 T = If,< 

T F ( A < T * I I ) .sci.c.li-d To 1 
r = A ( r . 1 1 ) 
no 2 1=1 *L 

2 a ( T • it =A ( I • i) /C 
1 rONTrwUE 

IF ( I T .FQ.Mr-O TO ^ 

1 . 1=1 * *1 
riQ 1 • T = J J * * 

IF< A f T * I I) . = 0.0. >t*0 To l'i 
no 3 '=1 *L 

o A(I*.n=A(II..O-A(I,Jl 

10 CONTINUE 

4 CONTINUE 
no 5 1=1 *K 

5 => I 1 ) - l> . 

P(K)-a(K*K+])/A(K.K) 

K K = K — 1 

no 6 1=1 *kk 
t=k-. i 

B(I)=A(I«K+1 ) 

LL = I + 1 
no 6 L =LL * K 

E> R(I)-H(I>-PfL>*A(I,L) 

RETURN 

FNU 
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SUHRniiriNE aINTFR 

COMMom ThETs ( 181 ) «U ( 1 <*1 > . I I , PI ,00.00 (ft'+O > •Mt7,iR).R(in,K 
nr mens i on r r 1 8i > 

00 = 0 . 

r '0 2 ■ J= 1 • 1 1 . 2 
ABERt!' . 
r < 1 > = 0 U > 

no. 3 T=2, J R1 
t hETr=ThH.T A» !>*»!/ I 80 . 

F < I) =0(11 *COS ( J*THETP> 

aRERrrA«EB+(F(l)+F(I-ll)/P. 

OQ ( J) =?.*ARFR/18f.i. 

return 

FND 


subroutine aonter 

COMMON THETA ( 181 ) ,Q ( 1°1 > , 1 1 .pi ,00 .00 (6S0> , A (RS.RM ,B (2S> ,K 
DIMENSION F ( 181 ) 

DO ? 1=1*11.? 

AREB=0 . 

F( 1 ); 0 ( 1 ) 
no 3 I =2 ,181 
THETp=THETA ( I>*PI/180. 

F(1)=0<I)*STNTJ*ThETR) 

AREB = aREBMF(I)+F(I-1i )/?. 

0Q< J)-2.*ARFR/180. 

return 

END 
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Input 

A single case consists of the determination of the pressure distribution around a 
given cone at a given Mach number and angle of attack. It is necessary to input the num- 
ber of terms in a Fourier expansion series to accurately represent the cross-sectional 
shape and the number of lifting and nonlifting solutions to be superimposed. For the 
loading routine used in the program, any column except the first may be used on the input 
cards. A decimal format is used for the input quantities. A description of the required 
inputs in the correct order and the FORTRAN variables used by the source program are 
as follows: 

Descrlption 


$NUM1 


MM 

N 

NN 

$ 

$NUM2 


M 

ALPHA 

$ 


Arbitrary name required by the loading routine to define the first input 
data block 

Number of terms necessary to accurately define cross-section geometry 
Number of nonlifting solutions to be superimposed 
Number of lifting solutions to be superimposed 
Denotes end of case (column 2) 

Arbitrary name required by the loading routine to define the second 
input block 

Free -stream Mach number 
Angle of attack, deg 
Denotes end of case (column 2) 


The system loading subroutine in the program (name list) is quite flexible in that 
the order of the input cards is unimportant and successive cases can be run by repeating 
the identification and $NUM cards followed by only the changed parameters and a $ card. 
An example of a set of input cards is given by the following listing of the inputs necessary 
to compute the pressure distribution around an elliptic cone (a = 0.57735 and b = 0.3) at 
M = 1.414 and a = 10°. 
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JjNU'-U 
K-1 
t j 

il>\l 

JO 

> i vl 0 *■'> ^ 


a L ^ H 

1> c! N 1 1 


^ . 

* J « 



Output 

An output listing for the example input is as follows: 


theta 

0 • 

1 .000000008+00 
2.00000000E+00 
3.00000000F+00 
4.00000000E+00 
5.00000000F+00 

6»oooooooof+oo 

7.00000000F+00 
8.000000008+00 
9.000000008+00 
l . oooooooqf+oi 
1 . 100000008+01 
1 . 200000008+01 
1 .300000008+01 
1.400000008+01 
1.500000008+01 
1.600000008+01 
1 .700000008+01 
1 .8U00000U8+01 
1 .90000000F +01 
2.000000008+01 
2.1000000GF+01 
2.20000000F+01 

2.3O0000u 0E + 01 
2.40000000F+01 
2.500000008+ul 
2.60000000F+01 
2. 700000005+01 
2.60000000F+Q1 
2.900000008+01 
3.00000000F+01 
3. 100000 OOF+Ol 
3.200000008+01 
3.30000000E+01 
3.4U000000F+01 
3.50000000F+01 

3.O0000000F+01 
3 . 7 0-0 0 0 0 0 0 E ♦ 0 1 
3.8U000000F+01 
3.9000000CF+01 
4.00000000F+01 
4. 1 00000 JOF+01 

4.20000000F+01 

4.30000000F+01 


CP* A=o 

5.10787096F-01 
5.10032787E-01 
5.07740580E-01 
5.038454088-01 
4.98288770E-01 
4.91068S29F-01 
4.82309418E-01 
4.72270014E-01 
4.613402328-01 
4.4949374GE-01 
4.38712Q43E-01 
4.279076728-01 
4. 178480 4 2E “01 
4.08629992E-31 
4.00182Q60E-01 
3.92317372F-01 
3.R4798509E-01 
R.77425984E-01 
O.70O96773E-O1 
7.82834780E-01 
R.5S7796H68-01 
7.491401 87F-01 
8. 43125946 E “01 
7,.37n79739£-01 
7 . 334897P8F-0 1 
3.29876410E-O1 
7.2641659OF-01 
7.23 19R7n3F“0 1 
7.20391 7?3E“0 1 
7. 17243766E-01 
7. 13924440E-01 
5 . 105209588-01 
7. 07206Q1 58-01 
7.04183490E-01 
7.01614771F-01 
7.99873433F-J 1 
2.980205578-Ul 
7.96ttl?7?4F-01 
9.95743783E-01 
7.94b0R?l IE-01 
2.932b6^05F-0 1 
9.91b3RO92F-01 
2.898156038-01 
9.879334Sf>E-01 


CP* A-CP* A = 0 

-2.7478366RE-01 
“2.741 1 04968-01 
-2.720906R4F-01 
-2.68729705F-0) 
-2.640555878-01 
-2.58) 403028-01 
-2.511 17846 E -01 
-2.431918978-01 
-2, 34626219 F- 01 
-2.25723695F-0I 
— 2. 16787628 F“ 01 
-2. 08O8364HF-01 
-1.9940S274E-01 
-1 .920536878— >1 
-1.84436471F-01 
-1 . 7R 0 8875 28-0 1 
-1 . 71+>90974F-o 1 
-1.655380R3F-01 
-l .596450696-01 
-1 .53SR50QOF-0 t 
-l .47O91772F-01 

-1 .4254bl62F-'H 
-1 .374497548-01 
-1 .3279165 IF-ol 
1 . 286 1 934 08-0 1 
-1 .249213028-01 

— 1 .21 k 2o296F — <• 1 
-1.185) 927 38 -01 
-1 .157693608-0 1 
-1.12 3619398-01 
-1 . 101 ?Sh6,)8-o] 
-1 . 07 ’4W3(i 28-0 1 
-l. 047733418-01 
-1.01 5849 1 38-01 
-9.897921698-02 
-Q c 66 3446198-02 
-9.458647948-02 
-9.28149R72F-02 
-9.125749178-02 
-8 . 9R v i 765 7 38-7? 
-R. R 74 1 R 4 4 7 8 — 0 ? 

-R. 687431 338-U? 

-8.516935908-0? 

— R. .34801 5 368 — 0? 


cptotal 

2.3600342R8-01 
2.3592J897E-01 
2.35649R96E-01 
2.351157048-01 
2.34231 1 R 38 — 0 1 
2.3292R227E-0 l 
2.311917728-01 
2.29j7R417F-U1 
2.2671401 78-01 
2.242700458-01 
2.21925314F-01 
2. 19R24024F-01 
2. 1H74276RF-01 
2. 165763098-0) 
2. 15 34644RF — 01 
2.142316208-01 
2. 1 31 u 75358-0 1 
2. 1 148790 18-01 
2.105517048-01 
2.091 496698-01 
2.07 7H7914F-01 
2.065940268-01 
2.056761 Q1 8-0 l 
2.05oRR0'7 6F-0l 
2.j4 w 1044-<8-01 
2. J4755104E-01 
9.o4 7402968 — 0 1 

2.0477944(18-01 

9.046218678-01 
2 • 04261 8268-0 1 

2. 0379858.08 — 01 
2.032726558-0 1 
2.U2R33574F-01 
2.02694557E-0I 
?.u26 351528-0) 
2 . 024384728-0 1 
9.034 340 7 7 E - 0 1 
2.034953418-01 
9. 0448629) 5— 01 

2.O460o6RR8-0 1 

?.O4'74«76O8-01 
2 . 0460 06 748-0 1 
2.O4616 04 68-T-01 
2.044735038-01 
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4.4UU00000E+O1 

4.5u0000008+01 

4.b0000000E+01 

4. fO000000F+01 
4.80000GG08+01 
4.90000000F+C1 
5.000000008+ni 
5. 10000000E+01 
5.2u000000F+01 
5.300000008+01 
5.4O000000F+01 
5.500000008+01 
5.600000008+01 

5. /uOOCOOOF+Ol 
5.80G000008+01 
5.90000000F + G 1 
6.00000000F+0) 
6 . 100000008+01 
6.200000008+01 
6.3000000CF+01 
6.400000008+01 
6.500000008+01 
6.60000000F+01 
6. 7U0000U0F+O1 
6 . 800000008+01 
6.900000008+01 
7.000000008+01 
7. 100000008+01 
7.2u0000008+01 
7.30000000F+01 
7 .h00000008+01 
7.SU0000008+01 

7 . bOOOOOOOF+O 1 
7 .700000008+01 
7.800000008+01 
7.900000008+01 
8.00000000F+01 
8 . 100000008+01 
8.200000008+01 
8.300000008+01 

8. “+00000008 + 01 
8.50000000E+01 
8 .b0000000F+0 1 
8.70000000F+01 
8.80000000F+01 
8.90000000E+01 
9.00000000F+01 
9.10000000E+01 
9.200000008+01 
9.30000000F+01 
9.40000000E+01 
9.50000000E+01 
9.60000000E+01 
9.70000000F+01 
9.80000000F + 01 
9. 90000000 E + 01 
1.000000008+02 
1 .01000000F + 0? 
1 .02000000E+02 
1. 0300000 0F + 02 
1.O4000000F+02 
1 .050000008+0? 


2.86181 1H7t-01 
?.84725a46fc-01 
7.836*37548-0 1 
7.82989*728-0 1 
7.825998658-01 
7.R?32?7068-01 
2.819703418-01 
7.813987828-01 
7.8oS456608-01 
7.794499298-01 
?.78?36*9bE-01 
7.77.,771 7^»E-01 
2.761 34R41E-01 
7.755147698-01 
2.75231FS6F-01 
7. 752^54408— Ul 

2 . 752874 8 3E—0 1 

2.753057428-01 

7.751 1 HO 348-01 

7.74t>54Qb5E~0 ) 

2.73b4079OE-01 

?. 73^819758-01 

7.72237b57E-01 

2.71566413E-C1 

?.7117bbl0F-01 

2.71O9S289E-01 

7.712588578-91 

P.71b33M0E-01 

7.717576058-01 

7.71790q2SF-01 

2.715573098-01 

2.71070849F-01 

•7.70426446E-01 

2.697738188-01 

7.692664258-01 

2.69017741E-01 

7.69061 834E-01 

2.693479468-01 

9.697520288-01 

2.701 17898E-01 

?.703054?68-01 

2.70234095E-01 

2.69907«98E-01 

7.69413874E-01 

7 . 6 8 894 7 8 7E -01 

2.685044898-01 

2.68359893E-01 

7.68504489E-01 

7.68H94787E-01 

2.69413874E-01 

2.69907898E-01 

?. 702340 95E“0 1 

2.703054268-01 

2.701 1 7898E“0 1 

2.69752028E-01 

2.69347946E-01 

7.69U61A34E-01 

7.690 1 7?4l E-0 1 

7.692664258-01 

7.69773818E-01 

7.70426*468-01 

7.71070R49E-01 


-ft. 16* 194? 18-0? 
-7.99424677F-0? 
-7.835929458-0? 
-7.69437355F-0? 
-7.S7-+28024F-O? 

-7. 47 a 4l954E- 0? 
-7.384555758-0? 
-7 . 30251 064F-0? 
-7.20*757f->2F-02 

-7.10779?6hF-0? 
-6.98 06071 8F-0? 
-6. 854847*48-0? 
-6.72 16428HF-0? 
-6.594441 74F-0? 
-6.48*5131 9F-<>? 
-6. 39] ftl587F-<V 
-6.31 1H5Q03E—)? 

-6.24786944F-0? 

-6.1 8 * 1 1 84648-0? 
-6.1 2 048 24.48- u 2 
— 6 .04416478F — ■)? 
-5 .95 u?? 0478-0 2 
-F .85?0?804F-.'j? 
—5 . 742878 1 98-02 
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Figure 12.- Comparison of results for elliptic cones without slender body assumptions. 

N = 32; M = \{2; a = 0°; a = 0.57735. 
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